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Abstract. In a recent paper [P. Mayer et al., Phys. Rev. Lett. 93, 115701 (2004)] it 
was shown, by means of experiments, theory and simulations, that coarsening systems 
display dynamic heterogeneity analogous to that of glass formers. Here, we present a 
detailed analysis of dynamic heterogeneities in the Glauber-Ising chain. We discuss how 
dynamic heterogeneity in Ising systems must be measured through connected multi- 
point correlation functions. We show that in the coarsening regime of the Ising chain 
these multi-point functions reveal the growth of spatial correlations in the dynamics, 
beyond what can be inferred from standard two-point correlations. They have non- 
trivial scaling properties, which we interpret in terms of the diffusion-annihilation 
dynamics of domain walls. In the equilibrium dynamics of the Ising chain, on the other 
hand, connected multi-point functions vanish exactly and dynamic heterogeneity is not 
observed. Our results highlight the similarities between coarsening systems and glass 
formers. 
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Introduction 

An obvious question to ask about glass-forming liquids is whether the dramatic slow- 
down of the dynamics on cooling is correlated with a corresponding increase in an 
appropriately defined length scale. Critical slowing-down around second-order phase 
transitions, for example, is correlated with the divergence of a static correlation length. 
In supercooled liquids and glasses, the consensus is that there is no growing static 
length scale, since the static structure — as measured, e.g., by the amplitude of density 
fluctuations — changes only negligibly while relaxation time scales grow by orders of 
magnitude. Any growing length scale in glassy systems must therefore reflect the 
spatial structure of the dynamics. In order for this spatial structure to be non-trivial, 
the dynamics must vary from point to point: it must be heterogeneous. The simplest 
conceptual picture of such dynamical heterogeneity is that some regions of a material 
are fast and others slow. The identities and locations of these regions may of course 
change over time. 

We will not try to review the literature on dynamical heterogeneity, which is 
vast, and refer instead to [U El 0]. We focus in this paper on the characterization 
of dynamical heterogeneities via multi-point correlations jU El El HI El EE CHI EH H2j ■ 
In the context of lattice models, one considers a spatial correlation function between 
sites % and j, of the general form = (FiFj) — (Fi)(Fj), but with Fi itself a two- 
time quantity such as Fi = Ai(t)Ai(t w ) and A some local observable. Thus CV,- is a 
four-point correlation function. It essentially measures how the relaxation of A at point 
i is correlated with the relaxation at point j. Such a definition should therefore pick 
up dynamical heterogeneities. Associated with Cjj is a so-called susceptibility, i.e., the 
spatially integrated correlation \ — ^ the length scale on which the 

dynamics is correlated grows, then so should this fourth-order susceptibility. 

It is worth noting that the four-point susceptibility x = can be re- 

expressed as x — (V-W) K<? 2 ) — (<?) 2 ] where q = ^2 i A i (t)A i (t w ) measures the "overlap" 
between configurations at time t w and t. From this we can easily gain qualitative 
insight into the dependence of x on At = t — t w . At At = 0, one does not expect large 
fluctuations in q; indeed, in a spin system and for A± = o~i, q is a constant for t = t w 
and hence x = 0- If? i n the opposite limit At — ► oo, the system decorrelates from its 
state at time t w , then one expects q to decay to a small value, again with only minor 
fluctuations. For intermediate times, however, the fact that a glass may remain trapped 
near the t w -configuration for a time that is both long and strongly dependent on the 
state at t w leads to large fluctuations of q between different dynamical histories and 
therefore to a large value of x- 

The usefulness of four-point correlations for understanding heterogeneous dynamics 
motivates us to consider their behaviour in coarsening systems jJJ. Because - contrary 
to glass formers - coarsening systems develop strong spatial correlations at long times, 
appropriately modified four-point correlations must be considered. It was shown in a 
recent letter that these reveal non-trivial spatio-temporal correlations in coarsening 
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systems which show strong analogies with aging glasses. This is, in fact, also the case for 
the Glauber-Ising chain, which we study in this paper. Its analytical tractability makes 
it an interesting candidate for exploring four-point correlations in coarsening dynamics. 
Results were partially announced in [T3|. 

Explicitly, the standard four-point function that is studied in the literature on 
dynamical heterogeneities becomes for our spin system 

Ci- k (At,t w ) = (o- k (t)a k (t w )o-i(t)o-i(t w )) - (o- k (t)a k (t w ))(ai(t)ai(t w )), (1) 

where t = At + t w > t w > and the times t, t w are measured from the quench. 
The definition Q implies that C n (At,t w ) is an even function of n. At At = it 
in fact vanishes for all n since Ising spins satisfy a 2 — 1. In the opposite limit 
At — > oo, because the configurations at t w and t will decorrelate, Ci~ k (At, t w ) approaches 
(o- k (t w )cri(t w )) (a k (t)ai(t)) - (o- k (t))(a k (t w ))(ai(t))(ai(t w )). The last term vanishes, so 
this limit reduces to the product of the spatial correlations at times t w and t — *■ oo. 
This argument holds quite generally in a spin system without an overall magnetization. 
However, in a typical glassy system the spatial correlations at times t w and t will be 
comparable and of limited range. In a coarsening system, on the other hand, the spatial 
correlations at time t have a diverging range as t — > oo, with {a k (t)ai(t)) — > 1, unless 
specific symmetries of the Hamiltonian are present |X2j|. This suggests that the large 
At limit of the standard four-point correlation C n (At, t w ) will be larger in a coarsening 
system than in glasses, with the growth of spatial two-point correlations obscuring 
genuine four-point correlation effects. 

This is indeed what we will find: to see genuine four-point correlations, we need to 
consider the "connected" four-point correlation 

Ci- k (At,t w ) = (o- k (t)a k (t w )ai(t)ai(t w )) - (a k (t)a k (t w )) (ai(t)ai(t w )) 

-(a fe (t)a i (t))(a fc (t w )a i (t w )) + (^(0^w))(^(*w)^(0>, (2) 

which differs from the standard version by the terms in the second line. One observes 
that the first of these just cancels the residual term from the four-point average in the 
limit At — > oo. We therefore expect that C„(At, t w ) — > for large At: the connected 
definition eliminates the uninteresting contributions from spatial two-point correlations. 
The second term in the second line of © can then be motivated as compensating for 
the first one at short times At, ensuring that, like C n , C n vanishes at At = 0. We note 
finally that C n is even in n as was the case for C n ; at n = 0, the definition (j2J) in fact 
implies that Co = at all times. 

We will also consider the four-point susceptibilities associated with JTJ), 0- These 
are defined as 

oo oo 

X (At,t w ) = J2 C n (At,t w ) and X(At,t w ) = ^ C„(At,t w ). (3) 

n=— oo n=— oo 

The layout of this paper is as follows. First we derive an exact expression for the 
connected four-point correlation C n (At, t w ) in Section [TJ In Section El the equilibrium 
behaviour of the standard and connected four-point functions is discussed. Scalings of 
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the connected four-point correlation and its associated susceptibility for non-equilibrium 
coarsening dynamics are analysed in Section^ We then interpret our results in Section^] 
in terms of the random-walk dynamics of domain walls. The behaviour of the standard 
four-point functions is briefly presented in Section 03 We conclude in the final Section. 



1. Derivation of the Connected Four-Point Correlation 



In this section we analyse the dynamics of the Glauber-Ising chain [T4] . quenched from 
a random initial configuration to some temperature T > 0. To recap briefly, the model 
has Hamiltonian 7i = —Y2i°'i <T i+i> where <7j = ±1 (i = 1, . . . N) are N Ising spins 
subject to periodic boundary conditions. Glauber dynamics consists in each spin oi 
flipping with rate \[\ — | 7 (7i(<7i_i + where 7 = tanh(2/T). 

General expressions for two-time multispin correlation functions in the Glauber- 
Ising chain are given in j!5j . for the finite model quenched at t — from equilibrium at 
an initial temperature 7] > to arbitrary T > 0. Let us now recall some results relevant 
for the present analysis: in the thermodynamic limit N — > 00 and for a quench from a 
random initial state T ; — > 00 we have the following representations 

(<7i(t)0i(t w )) = e- At I^At) +E$, 

{aiMvM°h{^)°h{^)) = [^Sf ] +H h - il (2At)]H ja . jl (2t w ) 



(4) 



-[ + e- 
+[-e- 



-At : 



4-^At) + £^}[ - e— 4_ i2 ( T At) + S s 



-At; 



■At : 



4_ i2 ( 7 At) + } ][ + , ( 7 At) + £ 



■ (ii,i2)i 
'12,32 J 



(5) 



for two and four-spin two-time correlation functions. In (jSJ) the indices must satisfy 
i\ < ii and j\ < 22- The general form of the functions £ and T is 

dimy) 

A %_ p ( 7 At)/4_ p (2t w ), (6) 



£ 3 



E 



T 3 



E 



p,q 



JJ sgn(j A - p) 

A=l 
dim(j) 

sgn(j A - p) sgn(j ; 

A=l 



g) 



e~ 2A * J, 



(7) 



In (JSj) and ((7|) the products over the sign-functions sgn(x), satisfying sgn(O) = and 
sgn(x) = o otherwise, involve the indices of all spins at the earlier time t w ; so when 
substituting (JUJ) into (J3J) there is only one factor sgn(j — p) while in (jSJ) we have two 
factors sgn(jx — p) sgn(j 2 — p)- The summations over p, q in (JBJ) , © are to be taken over 
the entire chain —00 < p,q < 00. Finally, the functions I q (t) denote modified Bessel 
functions ^H] while the H q (t) have the representation 

-t 



The physical meaning of if is 



9-1 



( 7 r) -7 9+ i( 7 r)] 



(8) 



(9) 
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but this holds for % < j only: in contrast to the two-spin correlation Q, the function 
H q (t) is odd in q and zero for q = 0. Further properties of H, which are summarized 
in are recalled below as and when required. 

Let us now focus on the connected two-time correlation Ci-k(At,t w ) defined in (J2J). 
In order to be able to express the four-spin term using (j^J) we require k < I; below, n 
always stands for I — k and is assumed to be positive. Also substituting (JU) and © for 
the corresponding two-spin correlations gives, after some rearranging, 

C n (At,t w ) = + H n {2At) - H n (2t)]H n (2t w ) 

-[<? - <f ] [<? + <f + 2 e" A * J ( 7 At)] 

+ +£gf ¥3 -£*S'° + 2e- A < I n (jAt)). (10) 

Here we have used £$ = £#> = £®, £$» = -E™ and E™ = -eff \ 

These properties follow directly from the definition of £ and reflect symmetries 
like (o~k(t)(Ti(t w )) = (cri(t)ak(t w )) . The problem of analysing C n (At,t w ) is now reduced 
to rewriting (|10|) in a convenient form. To this end one could utilize the closed 
representations for S and JF derived in [THj . These were, however, constructed for 
cases where the spins at the earlier time t w are at a fixed and small distance. In the 
current context this distance is given by n and we are interested in studying the scaling 
behaviour for n — ► oo or working out the infinite sum over n in Q. It is therefore 
necessary to develop a new approach for dealing with the expression (llOjl . 

As we show in the following it is convenient to rearrange the sums £ and T . Let us 
first consider the sums £ , appearing in (fTU|) only in very particular combinations. After 
substitution of (jOJ and a shift in the summation variable we obtain, for instance, 

<5 ± <f = E t 1 ± s ^ n + p)i s § n (p) e ~ At ( n ) 

This sum obviously reduces to a semi-infinite one due to the factor in the square brackets, 
for either choice of the sign. In order to lighten the notation and for the subsequent 
analysis it is convenient to introduce the weight function 

[ 1 q = 

w q = < 2 for 1 < q < n . (12) 
I 1 q = n 

We also use a slight modification of H 

H q (t) = 5 qj0 + H q (t), (13) 



and define 



S q = e~ At J Q ( 7 At) H q (2t w ) + 2 £ e~ A ' J p ( 7 At) H q+p (2t w ), (14) 

p=i 

oo 

T g = e" A * J n ( 7 At) ^(2t w ) + 2 e~ A * 4 +P ( 7 At) H q+P (2t w ). (15) 

p=i 
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In terms of (|12H15|) the two combinations of £'s in can then easily be shown to 
equal 

f( k ) _ C _ q- (-1 fl\ 

°k,k c k,k — ^n, l 1D J 

n 

4? + eSf + 2 e_A< J o(7At) = <S + $> 9 e~ A < 7,( 7 Af) ^(2t w ). (17) 

q=Q 

Similarly we find for the other two combinations of £'s in (I10J) 

n 

- <f } + 2e" A * 4( 7 At) = T + ^ ^ e" A ' J„_ g ( 7 At) # ? (2* w ). (19) 

Next consider the double sum T . Here the relevant combination is T^f* + H n (2 At) — 
H n (2t). Based on the identity I q (x + y) — J2 P ^p( x ) Iq+p(y) an d © one verifies that 

ff,_ fc (2t) -H t _ k (2At) = J2^' 2Atl k- P (lAt)I l ^( 1 At)H^ p (2t w ). (20) 
Expressing T via (|7J) and using (|2l]j) then yields 

T- fc ( J' + H n (2At) - H n {2t) = - [1 - s S n ( fc - P) s S n ( fc - 9) - p) s S n (^ " ?)] 

xe" 2A * 4_ p ( 7 At) I,_ 9 ( 7 A*) H q _ p (2t w ). (21) 

In analogy with (fTTj) the factor in the square bracket in ()21|) is non-zero only in a 
restricted range of the summation variables p, q. In fact, the two-dimensional infinite 
sum in (}2"Tj) may be rewritten as a finite number of one- dimensional semi-infinite sums. 
From this procedure, which is slightly cumbersome but trivial, and using the notation 
JH-HU) we obtain 

n n 

Ttf ] + H n (2At) - H n (2t) = e ~ At J </(7At) T ff - J> 9 e" A * J n _ g ( 7 At) S q . (22) 

q=0 q=0 

In terms of equations (jl6H19j) and (J2~2"j) our representation (fTUj) for the connected four- 
point correlation C n (At,t w ) is thus transformed into 

C n (At, t w ) — S n To — T n Sq 

n 

+ J2 w i e ~ At W At ) K H n (2t w ) - H q (2t w ) T n ] 

q=0 

n 

- J2 w o e ~ At In~ q {lAt)[S q H n {2t w ) - H q {2t w )S n ]. (23) 

q=0 

Equation (|2"3|) forms the basis for our subsequent analysis of C n (At,t w ). Note that up 
to this point we have not carried out any summations, except for (J2"U|) . The derivation 
of ()23)1 relies only on direct cancellations occurring in (|T0|) . 
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2. Equilibrium 

A striking feature of (|2*H|) is that it allows us to study the equilibrium behaviour of 
C n (At,t w ), still without working out any sums at all. To see this we first notice that in 
equilibrium we have, as can be shown ^H] from (JHJ) and (|T3|) . 

= lim H q {2t w ) = e for q > 0, (24) 

where 

i = ^ L = tanh 1/T . (25) 

7 

Due to the exponential dependence of H^ q on q, the expressions (|14|) . (J15)) for S, T 
satisfy in equilibrium 

<S^ = e<S eq and r; = fV. (26) 

It therefore follows immediately from (|2*H|) that our connected four-point correlation 
vanishes in equilibrium, i.e., 

C(Ai) = lim C n (At,t w ) = 0, (27) 

£ w — >oo 

for all n, At > and at any temperature T > as specified via 7. Looking back at the 
definition of C n (At, t w ) in equation (j2J), this implies an exact decomposition of four-spin 
two-time correlations into pairwise correlations. In other words, in equilibrium there are 
no genuine four-point correlations. 

Based on the general expressions given in ^2] we have verified that four-point 
correlations (o~k{t) cz^^w) o~i 2 (t w ) o~i 3 (t w )) likewise factorize. There should be a generic 
connection between this property and the fact that the Glauber-Ising model maps to free 
fermions J7j. However, as we will see below, the factorization only holds in equilibrium. 

Now consider for comparison the standard four-point correlation function It 
can be expressed in terms of the connected one, Equation (j2J), via 

Ci-k(At,t w ) = Ci^k(At,t w ) 

+ (a*(tM*)>MO*j(*w)> " (^fc(*)^(*w)X^(*w)^(*)>- (28) 
The first term vanishes in equilibrium, and so C^(At) = C n (At,t w — > 00) can be 
expressed purely in terms of two-point correlation functions. This suggests that the 
standard four-point function Ci-k{At,t w ) is strongly biased by pairwise correlations 
and thus not suitable for revealing genuine four-point correlation effects. 

3. Non-Equilibrium 

The significance of measuring connected four-point correlations becomes even clearer 
when considering non-equilibrium coarsening dynamics. For the sake of simplicity we 
focus on a zero-temperature quench, i.e., 7 = 1. The scaling of C n (At,t w ) in the limit 
of large times At, t w —>■ 00 and distances n — > 00 is then expected to be of the form 
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Formally, the long-time and long-distance limit is taken at fixed values of the scaling 
variables 

At \n\ , . 

a = - — and r] = — . (30) 



It 2*/r 

The first of these, a, measures the observation time interval At in units of the system's 
age t w while n is the ratio of distances n over the typical domain size, bearing in mind 
that the latter scales as O ( y^w) • The factors of | in the definitions of a, r\ are included 
for mathematical convenience in what follows. 

In order to obtain the scaling function fc(cn, rj) a leading order asymptotic expansion 
of (J23J) is required. To this end we use the asymptotic formula jTHj 

e^L(t)~-±=e-*'W, (31) 



/2%t 

which applies for q, t — > oo with q 2 /t fixed. In the same limit we have, by combining (JHj) 
with (}3*Tj) and the identity I q -i(t) — I q +i(t) = ?fl q {t), and setting 7 = 1, 

ff.«~»(4j) with iW4r*e-' (32) 



The function <&(x) is in fact just the complementary error function $(x) = erfc(x) = 
1 — erf (a;); we use the symbol $ to keep the notation compact. When substituting ([31)1 . 
(JH2J) into (JHJ), (fTHj) the sums defining iS 9 ,7^ turn into Riemann-sums such that 

2 ' ' " — » (33) 



(34) 





2v^w 

Using the expansions (|31 )) -(|34 )) in (J23|) and taking the scaling limit also turns the sums 
explicitly appearing in (J23*|) into integrals. In terms of our scaling variables a, rj the 
scaling function fc(a,n) is thus 

4 



f c (a,r}) = — \ du dve-^ u2+iv+v)2 } [$(r/ + - ®(u)<$>{r) + v)} 

.Jo Jo 

rr\ poo 

du dve-«l u2+{r > +v) \ [<S>(u + v)<S>(n) - ®(u)$(r) + v)} 
Jo Jo 

POO 

du / dve-^ v - u)2+v \ [$(u + v)$(rj) - $(u)$(r/ + v)) \ . (35) 



air 

+ 



Equation (J33|) is suitable for numerical evaluation; the resulting plots of fc( a ,v) are 
shown inFig. ^ One sees clearly that there are non-trivial four-point correlations in the 
non-equilibrium coarsening dynamics of the Glauber-Ising chain. We notice also that 
fc(ct, v) i s negative throughout (the plot shows the modulus); we will find an explanation 
for this feature in Section 0] below. Because of the rather rich structure in fc{a,n) let 
us next discuss its various scaling regimes. 

First consider a 1, which corresponds to At 3> t w . The behaviour of fc(oc,T]) in 
this regime is easily obtained form (|3*K|) by Taylor expanding the exponentials in I /a. 




Figure 1. Modulus of the scaling function fc{a,rj) of the connected four-point 
correlation C n (At, t w ) versus distance 2r\ = \n\/y/t^ for various time ratios 2a = At/t w . 
The curves are obtained by numerical evaluation of the exact scaling functions 
The dashed-dotted straight lines are the asymptotes (f^|> and where the sloping 
line corresponds to (|42[) for the regime ry 2 <gC a <C 1 while the horizontal ones represent 
lfi3|) with 2a — 10~ 3 , 10~ 2 and apply in the regime a < (j 2 < 1. The data is discussed 
below and interpreted in Sectional 



The zeroth order contributions, where the exponentials are replaced by unity, vanish: 
in the first integral in (f33|) the combination of $'s is antisymmetric under exchanging 
u,v while the second and third integral cancel. So the leading behaviour of fc(a,rj) 
follows from first-order contributions where, for the same reasons, various terms drop 
out. The remaining integrals, which are of the type J dxx l Q(x) with i = 0,1,2 and 
f dx f dy (x + y) $(x + y), are solvable [211 and lead to the large a expansion 

fc{a, V ) = --^ { 12 v 2 e"" 2 - 2 [(1 + Q V 2 ) - (1 + tjVI 0Ftj $(77) 

-(3 + 2?7W^)}^ + 0(^. (36) 

This means that the magnitude of four-point correlations drops as I /a 2 for a ^> 1, at 
all distances r\. The plots of the exact fc(a,rj) in Fig. [T]for 2a = 10 1 , 10 2 illustrate this 
nicely. Details of the shape of fc(&,v) a t large a are revealed by expanding (|3l)j) for 
i] ^ 1 and T] 3> 1, i.e., distances n far below and far above the typical domain-size at 
time t w . In the former case we simply Taylor expand in rj while in the latter we use ^E] 
that = 1/(^/717]) e~' q2 \l — \/(2r] 2 ) + 0(r]~ 4 )]. The leading terms are 

f c (a, rj) ~ -- (- - l) ^ for r, 2 < 1 « a, (37) 

7T \ 7T / CT 
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fc(a,r))~--^e-*> 2 for l«a, 77 2 . (38) 

Again the plots in Fig. Qfor let = 10 1 , 10 2 clearly show the power-law behaviour given 
by (|S7jl and the Gaussian cutoff (|5Sjl. 

Now we turn to a <C 1 or At <C t w . To obtain expansions of fc{&, v) i n this regime 
it is necessary to rearrange (|35j). It is further convenient to introduce 

fc(<x,p) = f c (a,y/ap) with p = — = — _ . (39) 

y/a V2At 

The scaling variable p may be viewed as an alternative measure for distance and thus 
replaces rj. In this notation and by shifting and scaling the integration variables equation 
(l33|) can be re-expressed as 
4 f f°° f°° 

/c(a,p) = -|y Dm y [$ a (p + u)$ Q (t; -p) - $ a (u)$ a (i;)] 

+ Du Dv [$ a (v + u- p)$ a (p) - ® a (u)<f> a (v)} 
Jo J p 

Bu J Dv - it + p)$ a (p) - $ a (p - n)$ a (p + v)] j . (40) 

Here we have introduced the short-hands Dx = e~ x2 dx and <& a (x) = $ (y/ax). Equation 
(jlQI is suitable for studying the small a regime of fc(a,rj). Two cases have to be 
distinguished: we can expand around a = either at fixed p or at fixed rj. When 
keeping p fixed we effectively look at distances rj = p y/a = O(yfa), while obviously 
i] = 0(1) if we fix 7]. Because the p and r\ length scales become disparate for a — > 0, 
separate expansions must be made. 

The shape of fc(a,n) for small a and fixed p immediately follows from (f4T)j) by 
Taylor expanding the functions <& a (x) = <&(y/ax) in y/a. This turns the integrands 
in (jHH) into Gaussians (contained in Dm, Dv) with polynomial factors. Evaluating the 
integrals gives 

/ c (a,p) = --i{i_ [e^ 2 + v^p(2-$(p))] [e^ 2 -v^p<l>(p)]}a + 0(« 3 / 2 ). (41) 

To understand the result (f4*T]) it is instructive to consider the limits p <C 1 and p 3> 1, 
corresponding to distances n <C y/At and n ^> y/At, respectively. One finds to leading 
order and in terms of a and n 

fc(a,r ] )~--(l--)r ] 2 for n 2 < a < 1, (42) 



7T \ 7T 



4 

/c(a, 77) ~ j a for a«?| 2 « 1. (43) 



7T 2 



We note that at small p the leading term in (|41j) is p 2 a = rf and thus (|42|) follows. So 
fc{&, rj) initially grows as rj 2 independently of a. But for rj w a/ck four-point correlations 
level off at a plateau of height 0(a), as given by equation (JJ3J). The asymptotes (J4*2"j) . 
(|43|) are shown in Fig. ^ 

It remains to discuss the behaviour of fc(c*,v) f° r small a and fixed n. In this 
case p = 77/ y/a in (|3U|) diverges for a — > 0. Because of the Gaussian weights in Dm, Dv 
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only integrals containing the neighborhood of u = v = then contribute significantly to 
fc{ct,rj). This holds for the third integral in (|40p. but not for the other two. The first 
one, for instance, satisfies the bound 
4 

7T 



POO poo 

Vu Dv [$ a (p + u)$ a (v-p) -$ a {u)$ a {v)] 

JO J p 



< 2$(77)$(p). 



This follows from the triangular inequality, the identities = f°° Dx and 

$(77) = $ a (p) and the fact that is monotonously decreasing. The same bound 

can be used for the second integral in ()40|) . Extending the w-integration range in the 
third integral to J °° Dm likewise only produces excess contributions of the same size. 
Therefore, to order O ($(^)$(p)) equation (j4*Uj) reduces to 

j POO POO 

fc{a, 77) ~ / Du Dv [$ {^fa{v - u) + 77) $ (77) - $ (77 - y/au) $ (77 + 

71 Jo Jo 

at small a and fixed 77. Here we have substituted = ^(y^x) and yfa p = 77. 

As will become clear in a moment the above expression has power-law scaling at small 
a. On the other hand, $(p) ~ \fa/{\fTi 77) e™" 2 /" vanishes faster than any power-law 
for a — > 0. Therefore we may safely ignore the O ($(7/)$(p)) contributions discarded 
above. By Taylor expanding the last representation for fc(ct,v) i n V®i which leads to 
simple Gaussian integrals, we thus finally obtain the small a scaling at fixed 77, 

f c (a, 77) = --L^ 2 fe-" 2 - V5Ftj $(77)! a + O (a 3/2 ) . (44) 

Equation (jUj) tells us that, on the length scale set by 77, four-point correlations decrease 
linearly with a as a — > 0, for any value of 77. The plots for 2a = 10~ 2 , 10 -3 in Fig. ^ 
illustrate this. At small and large 77 = |n|/(2y / t^) the behaviour of flU)) is to leading 
order 

4 

f c (a,r})~ -a for a<7/ 2 <l, (45) 

/ c (a,77) ^e- 2 ^ 2 for a « 1 « V 2 . (46) 

TV Tj 

The plateau obtained at small 77 matches the one at large p, Equation (|35Jl . as it 
should. For small a and large 77, finally, /c(a, 77) has an essentially Gaussian cutoff 
The latter occurs at slightly smaller 77 than in the large a case, see ()38|) and Fig. ^ 

Our discussion of /c(a, 77) has given us a complete understanding of the spatio- 
temporal scaling of the connected four-point correlation C n (At, t w ). The various scaling 
regimes in Fig. □ are characterized by JED, (JSHJ) , (E2), (ESI), (US) and (gEJ. We 
now conclude our analysis of the connected four-point correlation in nonequilibrium 
coarsening by studying the scaling of the associated four-point susceptibility. From its 
definition Q, the properties C_ n (At,t w ) = C n (At,t w ), C (At,t w ) = and the scaling 
(|29|) one finds in the large-time limit At, t w — > 00 at fixed a, 

00 00 (At n \ f°° 

^(AMw) = 2VC n (AMw)~2V/ c — ,— = ~4vC / d77/ c (a,r7). 

n=l n=l \ Mw Z VW Jo 
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This defines the scaling function Fcipt) for the susceptibility, 

X{te,U)~y/UFc(^f\ with F c (a)=4jf dr)fc(a,r)). (47) 

Numerical integration of (J4T|) using fc(oc, rj) as given in produces the plots in Fig.EJ 
Because fc(a,r)) < we also find that F c (a) is negative throughout. Furthermore, 
Fc{o) shows power- law behavior at small as well as large a: 

8 (\/2 — l) 

Fc{oc) = , /9 ) a + 0(a 3/2 ) for a < 1, (48) 

, 2 (8^-9) 1 / 1 \ 

Fc(a) = - > -, + - for a » 1. (49) 



57T 3 / 2 a 2 V^ 3 , 
The expansion for large a is obtained by substituting the corresponding expansion (J3fij) 
of fc(a,rj) in (|47|). At small a we split / °° d^ = dr] + J^°^dr] in (J4T|) with some 
^ > 0. The integrals correspond to the p and rj length scales where the expansions 
(}4Tj) and (|44p of fa(a,r]) apply, respectively. One easily shows that the contributions 
to -Fc(qO from the p length scale are 0(a 3 ^ 2 ) while those from the r\ scale grow as 
0(a). Therefore the leading term in the small a expansion (|48jl is given by substituting 
()44|) into (|47j) . The fact that -^(tt) vanishes at small and large a is not surprising: 
as discussed below (J2J) the connected four-point correlation C n (At,t w ) goes to zero for 
At — ► as well as At — > oo. However, that the approach to these limits is through the 
power laws (}4~%j) and (J4T?|) is rather non-trivial. From the arguments above it is clear 
that the linear scaling at small a is a consequence of the same scaling of the plateau 
height in fc(oc,r]). We consider next a physical picture which provides some intuition 
into how this scaling arises. 



4. Random Walk Interpretation 

The exact scaling results summarized in Figure have a lot of structure. To develop a 
physical understanding, we now use the fact that low-temperature dynamics of the Ising 
chain can be described in terms of domain walls that perform independent random walks 
with diffusion rate | until they meet, when they annihilate with rate close to unity |18j . 
In equilibrium there is also the reverse process, where pairs of domain walls or "walkers" 
are created at a small rate. 

In a space-time diagram jS], Figure E^a), the spins <7fc(t w ), 07 (£ w ), a k {t) and &i{t) 
define the four corners of a rectangle. Spin products are then determined by whether an 
even or odd number of domain walls cross the relevant edge of the rectangle. For instance 
o~k{tvi)o~k{t) equals 1 if an even number of walkers cross the bottom edge; otherwise it 
equals —1. One can then classify all possible situations by the parity of the number 
of walkers crossing the four edges. Numbering the edges in the order left - right - top 
- bottom, we will for example denote the situation where an odd number of walkers 
crosses on the left and bottom as 1001. Because walkers only annihilate or recreate in 
pairs, the total number of walkers crossing the rectangle has to be even, so that there 
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Figure 2. Modulus of the scaling function Fc(a) of the four-point susceptibility 
X{At,t w ) defined in gJJ, versus 2a = At/t w for zero temperature coarsening. Dashed 
lines represent the asymptotes l|48[) and (|49|l . Inset: Fc on a linear scale. 



are eight possible situations. After a short calculation one shows that, in terms of the 
corresponding probabilities, the connected four-point correlation is 

C n (At,t w ) = 8(pooooPiin +P0011P1100 -P0101P1010 -PonoPiooi)- (50) 

The last two terms are always identical to each other, due to the spatial mirror symmetry 
of the problem. The representation (fKUj) is particularly useful when the number of 
walkers crossing the rectangle is small. We expect this to be the case in the regime 
n l/c(t w ) and y/At l/c(t w ), where c(t w ) is the concentration of domain walls at 
time t w . In the coarsening regime, where c(t w ) ~ 1/ (2y/ir t w ), these conditions translate 
to n 2 t w and At <^t w . In equilibrium, on the other hand, c is a constant fixed by the 
temperature. To include both cases, we simply write c below, even for the coarsening 
case. 

Our restrictions so far still leave open the relative magnitude of At and n 2 ; let us 
focus first on the case where At <C n 2 . The space-time rectangle is then wide in the 
space-direction and narrow in the time-direction. From this one can deduce the leading 
contributions to the various probabilities, which are shown in Figure 0(b). We have 
denoted by T the probability that a walker will cross a corner of the rectangle. Since 
n ^> \/At, this probability is dominated by the smallness of y/At, while the spatial 
extent of the box is irrelevant. Thus, T can be calculated as the probability that a 
walker will cross from one half-space into the other during the time interval At, which 
is easily found as T = c\J At/2ir. The other quantity that appears in the probabilities is 
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Figure 3. Space-time representation of the connected four-point correlation; grey 
lines indicate domain- wall trajectories Panel (a): snapshot of trajectories over a 
spatial region (vertical) of 500 sites and a time-window (horizontal) of At = 1500 for 
equilibrium dynamics at T = 0.6. Panel (b), (c): schematic trajectories. See text for 
discussion. 



the factor r, which we define to be the joint probability at time t w of having two walkers 
at distance n, normalized by c 2 . With this normalization, one would have r = 1 for 
n ^> 1/c because correlations between walkers vanish at large distance. At equilibrium, 
where walkers are uncorrelated at any distance, one in fact has r = 1 even for n < 1/c. 
In the coarsening situation, on the other hand, r is of order nc for n C 1/c and thus 
vanishes to leading order. This reflects the effective repulsion between walkers: a walker 
that has survived the coarsening dynamics up to t w is not likely to have other walkers 
within a distance of order 0(y/t^) = 0(l/c). Putting the terms from Figure Efb) 
together gives 

C n (At, t w ) « 8(1 x 2PV + 2PV x nc - 2T 2 ) w 16r 2 (r - 1) = -(r - l)c 2 At, (51) 

7T 

where we have used the fact that nc < 1 to neglect the second term. This simple 
expression explains two important qualitative observations made above. First, in 
equilibrium, C n vanishes because r — 1. Second, in the coarsening case, C n is negative 
because r < 1 when n 1/c. Thus, the sign of C n arises from the effective repulsion 
between walkers discussed previously. Quantitatively, using that r < 1 for n C 1/c, 
the result (J5Tj) predicts that 

8 9a 2 At 4 
C n «— c 2 At - — = -a. 52 

7T 7l Z t w 7l Z 

This is precisely our expansion (f4*Hj) for the regime At < n 2 < where C n has an 
rz-independent plateau whose height increases linearly with At/t w . So we now have a 
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microscopic picture for the occurrence of this plateau in terms of domain-wall dynamics. 

Next we apply similar arguments to the regime where n is small compared to y/At, 
n 2 -C At < t w . As shown in FigureEfc), the space-time rectangle is now extended in the 
time-direction. As a consequence, the two probabilities poon and pnoo swap their leading 
contributions. For poon, the leading term is now produced by a single walker crossing 
the rectangle from top to bottom or from bottom to top. Since n <C \/At, the spatial 
width of the box can be neglected to leading order, and Poon reduces to the probability 
of crossing from one halfspace to the other. Bearing in mind that the crossing can occur 
from the top or the bottom then gives poon — 2c\/At/2~K. For pnoo, one might naively 
expect to get a product of two corner-crossing probabilities times a repulsion factor. 
However, as opposed to the case of pnn, the two corner crossings shown can in fact be 
achieved by a single walker which starts and ends within the interval of size n. This gives 
a leading contribution of pnoo = nc x n/V^nAt. The remaining terms are calculated 
as before, except for the fact that now the corner-crossing probability is T = nc/2. 
Assembling all terms, and using again that r 1 to neglect PooooPim = 2T 2 r, we thus 
get for the coarsening case 



C n (At,t w ) ^8 




At n 2 c n 
2tt v^At 



. (53) 



This again has the correct negative sign overall. It also predicts that, in this small- n 
regime, C n grows quadratically with n, with an amplitude independent of At. In fact, 
using n 2 c 2 ~ n 2 /(47rt w ) = rf/n, the result (JHSJ) coincides with the expansion (jl^jl as it 
should. 

The random walk picture has turned out to be useful for explaining the behaviour 
of C n when n 2 , At -C t w . In the remaining regimes discussed in the previous Section, 
on the other hand, where either n 2 or At are large compared to t w , it is less helpful 
because a large number of annihilating walkers has to be considered. It is then no longer 
obvious how to estimate the probabilities in ()50|). Nevertheless, the Gaussian cutoff for 
n 2 ^> t w that we found in (jSHJ), is at least qualitatively reasonable: for length scales 
n ^> \/tZ, correlations between random walkers become weak and one should effectively 
retrieve the equilibrium situation, where C n = 0. 



5. Standard Functions out of Equilibrium 

In Section |21 we saw that while the connected four-point correlation C n (At,t w ) and 
its associated four-point susceptibility ^(At,^) vanish in equilibrium, the standard 
functions C n (At, t w ) and x(At, t w ) are biased by two-spin correlations. It is the 
purpose of this section to show that the same is true for the non-equilibrium coarsening 
dynamics. The link ()28|) between the standard and connected four-point correlations 
allows us to express their difference AC„(At,t w ) = C„(At,t w ) — C„(At, t w ) purely in 
terms of two-spin correlations. This makes the analysis of AC„,(At,t w ) rather simple: 
spatial correlations are given in (jHJ) in terms of H n while temporal correlations for zero 
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temperature coarsening have the exact representation |T3] 



2f„ 



M*)<ri(*w)> = e" (i+iw) ^ 4(* + t w )+ drl n (t + t w -r) [J (r) + h(r)) \ . (54) 



In the scaling limit At,t w ,n — > oo with a = At/(2i w ) and 77 = | rz | / (2 v^w) fixed we 
substitute the expansion (J3TJ) into ()54p. Combining terms according to (|2*8Jl and some 
rearranging then produces AC n (At,t w ) ~ f^ct^ri) with 

/ ,. N I /•arccoty'a ^ 2 

f A (a, rf) = $(77) $ - 2 I - ^ - / dz e~^ 1+ ^ scc2 W I , (55) 



where sec(2;) = 1/ cos(z). The scaling of the difference between the four-point 
susceptibilities A%(At,t w ) = x(At,t w ) — X(At, t w ) = AC n (A£, i w ) then follows 
by analogy with (|37j) : writing A%(At, t w ) ~ v / ^"-Pa(«) we obtain via integration of 
(J55J) over 7], 



F A (a) ~ -4= < VI + 2« - 1 + - 

% 7T 



/ — / 1 + a I — J r 

arctan V2a — \ arctan 2 y «(1 + a J 



(56) 

Plots of f&.(a,,rj) and -Fa(«) are shown in Figure 0J Comparing the vertical scales in 
Figures HI IH demonstrates that the standard functions C n (At,t w ) and x(A£, i w ) are 
completely dominated by the two-spin contributions (|55j). (|56)l. A plot of the four- 
point susceptibility Xn(At,t w ) = X n (At,t w ) + Ax n (At,£ w ) ~ V^[-^c(a) + ^a(o;)], for 
instance, would be indistinguishable by eye from the inset of Figure EJ 

Therefore, as claimed, the standard four-point function (0) and its associated four- 
point susceptibility (j3J) are not suitable for measuring genuine four-point correlations 
in the coarsening dynamics of the Glauber-Ising chain. In comparison to strongly 
heterogeneous systems the relative magnitudes of, e.g., the connected four-point 
susceptibility (inset in FigureEJ) and the corresponding two-point bias (inset in FigureEJ) 
are reversed in coarsening systems. 



6. Conclusions 



In this paper we have explored dynamical heterogeneities in coarsening systems by 
studying multi-point correlations in the dynamics of the Glauber-Ising chain. Since 
conventional four-point correlation functions become dominated by strong spatial 
correlations that develop in coarsening systems at late times, we considered "connected" 
four-point functions where these uninteresting two-point contributions are eliminated. 
We were able to obtain exact results and scaling forms for these functions and the 
associated spatial integral, i.e. the connected four-point susceptibility. As a function of 
the time difference At, this multi-point susceptibility has an extremum, as is found in 
glass formers, for times of the order of the waiting time, indicating the timescale for 
which dynamic heterogeneity is maximal. 

Interestingly, we found that the connected four-point susceptibility is negative 
throughout, and we were able to give an interpretation for this behaviour in terms 




Figure 4. Plots of the scaling expansion i|55|) of AC n (At,t w ) versus scaled distance 
2jj = \n\/y/t^ for various time ratios 2a = At/t w . Inset: Normalized contribution 
F/\(a) of the two-spin terms to the four-point susceptibility x(At,t w ), Equation l|5t)[) . 



of the dynamics of domain walls, which undergo free diffusion and pair annihilation. 
The negative sign of the susceptibility directly reflects the fact that there is an effective 
repulsion between the domain walls, each having a "depleted zone" around it where 
the likelihood of finding another domain wall is low. At equilibrium, on the other 
hand, domain wall positions are uncorrelated and this leads to the vanishing of the 
susceptibility, and of the underlying four-point correlations, for all At. This latter result, 
which we established using explicit expressions for four-spin correlations, appears rather 
non-trivial. It would be interesting to verify whether it also extends to equilibrium 
correlation functions of higher order. If it does, one suspects that there should be a 
deeper reason, possibly related to the mapping of the Glauber-Ising chain dynamics to 
free fermions |17j . 

We also discussed the spatial dependence of the connected four-point correlation 
functions. This has a richer structure than one might have expected, but the random 
walk picture again gave a good qualitative (and, in some regimes, quantitative) 
understanding of our exact results. In future work, it will be interesting to see if and 
how our findings generalize to other coarsening systems. Encouragingly, simulations 
show that many of the key features we found here extend at least to two-dimensional 
Ising models [T5] . 
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